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Quantum transport in inhomogeneous magnetic fields is investigated numerically in two- 
dimensional systems using the equation of motion method. In particular, the diffusion of electrons 
in random magnetic fields in the presence of additional weak uniform magnetic fields is examined. 
It is found that the conductivity is strongly suppressed by the additional uniform magnetic field 
and saturates when the uniform magnetic field becomes on the order of the fluctuation of the ran- 
dom magnetic field. The value of the conductivity at this saturation is found to be insensitive 
to the magnitude of the fluctuation of the random field. The effect of random potential on the 
magnetoconductance is also discussed. 



I. INTRODUCTION 



Quantum transport in two-dimensional (2D) disor- 
dered systems in inhomogeneous magnetic fields has been 
studied extensively, in connection with the fractional 
quantum Hall effect (FQHE)J fn particular, the An- 
derson localization 2 in a 2D system in random mag- 
netic fields (RMF) with zero mean, which arises in the 
mean field theory of the fractional quantum Hall ef- 
fect at filling factor v — 1/2, has been studied by 
many authors From the theoretical point 
of view, whether such a system has a metallic phase or 
not has been an important issue. Systems in random 
magnetic fields belong to the 2D unitary universality 
class which is, in general, expected to have no metal- 
lic phase. If the system has a metallic phase, it throws 
a question on the validity of the conventional classifica- 
tion of the universality classes of the Anderson transition. 
Although many numerical as well as analytical attempts 
has been made to clarify this point, the conclusions are 
so far still controversial. On the other hand, the singu- 
lar behavior of the conductance fluctuations^^ and the 
density of states at the band center (E — 0) has been 
investigated and argued that these singularities are gov- 
erned by the chiral unitary universality class ^*i4 

Apart from such a theoretical interest, the transport 
properties of 2D electron systems in random magnetic 
fields are important to understand the experiments in 
the quantum Hall systems. In fact, the structure ob- 
served in the magnetoresistance at the filling factor 1/2 
has been analyzed using models with random magnetic 
fields^ On the other hand, recently, two-dimensional 
electron systems in random magnetic fields have been 
realized experimentally by attaching small magnets on 
the layer parallel to the 2D electron gas in a semicon- 
ductor heterostructureii^ In such a system, one can also 
control the strength of the random magnetic fields. The 
magnetotransport has been measured and interestingly, 
it is found that the magnetoresistance in systems hav- 
ing random magnetic fields exhibits similar structures to 



that in the fractional quantum Hall system at the filling 
factor v = l/2>i& Since the origins of the magnetic field 
fluctuations in these two systems are quite different, it 
is important to clarify the origin of this similarity in the 
magnetoresistance between these two systems. 

One of the significant features of the magnetoresis- 
tance observed in the recent experiment^ is the increase 
of the magnetoresistance in two different scales of the 
magnetic field. One is the dip structure around zero 
field whose width is on the order of the fluctuation of 
the random field, and the other is the positive magne- 
toresistance in a larger scale of the magnetic field fol- 
lowed by the Shubnikov-de Haas oscillations. In recent 
years, the semi-classical theory for the magnetotransport 
in a smoothly varying random magnetic field has been 
developed^ where the snake state near zero magnetic 
field lines plays an important role. In the case of a weak 
disorder, a pronounced positive magnetoresistance com- 
ing from a classical origin has been reported for small 
magnetic fieldsJ^ Although the semi-classical approach 
provides qualitatively consistent results with the exper- 
iments, it is also important to examine the transport 
property by the quantum mechanical calculations for its 
further understanding. 

In the present paper, we study the magnetoconduc- 
tance in the random magnetic field in two dimensions. 
We consider a tight-binding model in random fluxes hav- 
ing no spatial correlation, which has been analyzed by 
many authors for the study of the Anderson localization 
in random magnetic fields. In two dimensions, it has 
been shown that such a model has the large localization 
length near the band center and the diffusive behavior is 
observed in that energy regime for numerically accessi- 
ble length scales. 11 Since there is no spatial correlation in 
random magnetic fields, the random magnetic field is not 
smooth and hence the validity of the semi-classical con- 
cepts is not obvious for the present model. In our model, 
we find that the conductivity is strongly suppressed by 
adding a uniform field and also that it shows a satura- 
tion when the uniform field becomes on the order of the 
fluctuation of the random magnetic field. 
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We adopt the equation of motion method in order to 
examine numerically the diffusion of an electron in in- 
homogeneous magnetic fields. Adopting this method, we 
have an advantage that very large systems compared with 
other numerical methods can be considered. However, 
within the present method, we are able to estimate the 
longitudinal conductivity (<j xx ) only, and not the Hall 
conductivity {c xy ). Due to this, we are not able to discuss 
the magnetoresistance directly. As mentioned above, the 
2D random magnetic field system has been shown to have 
very large localization length near the band center. It is 
then natural to assume that, in the present model, this 
regime is responsible for the metallic behavior observed in 
experiments. We therefore confine ourselves to the case 
that the Fermi energy lies near the band center. First, we 
discuss the transport in the random magnetic fields with 
zero mean and next, we examine the effect of additional 
uniform magnetic fields. The strength of the additional 
uniform magnetic field considered in the present study is 
assumed to be on the same order as or smaller than that 
of the random magnetic field. Implications from our nu- 
merical results for the magnetoconductance are discussed 
in comparison with the recent experimental results. 



II. MODEL AND METHOD 

In order to describe the two-dimensional system in ran- 
dom magnetic fields, we consider the following Hamilto- 
nian 

H = VeM^)C}C 3 + £ £<C?C< (1) 

<i,j> i 

on the square lattice. Here Cj(Cj) denotes the cre- 
ation(annihilation) operator of an electron on the site 
i and {e^} denote the random potential distributed inde- 
pendently in the range [—W/2,W/2]. The phases 
are related to the magnetic fluxes {<j>i\ through the pla- 
quette (i, i + x, i + x + y, i + y) as 

6i,i+x+6i+x,i+x+y +6i+x+y,i+y-\~9i+y,i = — 2,TT(f>i/<j)o (2) 

where (f>o — h/\e\ stands for the unit flux. The fluxes 
{4>i} are also assumed to be distributed independently in 
each plaquette. The probability distribution P(<f>) of the 
flux cj) is given by 

The variance of the distribution is accordingly given by 

(Mi) = %<t>l^, r (4) 

In order to solve the time-dependent Schrodinger equa- 
tion numerically, we employ the decomposition formula 



for exponential operators— The basic formula used in 
the present paper is the forth order formula 

exp(x[A 1 +- • -+A n ]) = S(xp) 2 S(x(l-4p))S(xp) 2 +0(x 5 ), 

(5) 

where 

S(x) = e xAl/2 ■ ■ • e^-i/V^e^- 1 / 2 • • • e xAl/2 . (6) 

The parameter p is given by p — (4 — 4 1 / 3 ) -1 and 
A\,..., A n are arbitrary operators. We divide the Hamil- 
tonian into five parts as in the previous paper— so that 
each part is represented as the direct product of 2 x 2 
matrices. By applying this formula to the time evolution 
operator U(t) = exp(— iHt/U), we obtain 

U(St) = U 2 (-i P 5t/h) 2 U 2 (-i(l-4p)5t/fi)U 2 (-ip5t/H : f) 
+0{6t 5 ) (8) 

with 

U 2 (x) = e xIi ^ 2 • .. e *"VV ff5 e xi W 2 • • -e xH ^ 2 , (9) 

where H = Hi H + H$. It is to be noted that the U2 

can be expressed in an analytical form while the original 
evolution operator U can not be evaluated exactly with- 
out performing the exact diagonalization of the whole 
system. We are thus able to consider larger system-sizes 
than other numerical methods, such as exact diagonal- 
ization and the recursion method based on the Landauer 
formula, which is one of the advantages of the present 
method. This method has already been successfully ap- 
plied to the case of W = and h T i = 1— as well as to 
the 2D symplectic class— 

The system we consider is the square lattice of the size 
999 x 999 with the fixed boundary condition. All the 
length scales are measured in units of the lattice con- 
stant a. To prepare the initial wave packet with energy 
E, we numerically diagonalize the subsystem (21 x 21) 
located at the center of the whole system and take the 
eigenstate whose eigenvalue is the closest to E as the ini- 
tial wave function. In the following, we set the energy 
E to be E/V = —0.5 which is close to the band cen- 
ter. The single time step 5t is set to be St — 0.02^/^. 
With this time step, the fluctuations of the expectation 
value of the Hamiltonian is safely neglected throughout 
the present simulation (t < 200^/^). We observe the 
second moment defined by 

(r 2 (t)) c ^(r 2 (t))-{r(t)) 2 (10) 

with 

(r"(0)=5>»bKr,t)| 2 , (n=l,2) (11) 

r 

where ^>(r, t) denotes the wave function at time t. In the 
diffusive regime, the second moment is expected to grow 
in proportion to t 

(r 2 ) c = 2dDt, (12) 
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FIG. 1: The second moment as a function of time for ft r f = 
0.2, 0.25, 0.3, and 0.4. Five realizations of random magnetic 
fields are considered for each h T f. 
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FIG. 2: Conductivity in units of e 2 /h for h lf = 0.2, 0.25, 0.3 
and 0.4. 



where the diffusion coefficient is denoted by D and d 
is the dimensionality of the system. The diffusion coef- 
ficient D is related to the conductivity by the Einstein 
relation a = e 2 Dp. 2Q Here p denotes the density of states. 
In Fig. 1, the second moment is shown as a function of 
time for the case of W = 0. For each value of h T {, five 
realizations of random magnetic fields are considered. It 
is clearly seen that for t > 50h/V the second moment 
increases in proportion to t, which means that the sys- 
tem is in the diffusive regime. For t < 50h/V, we see 
the ballistic behavior, where the second moment grows 
as oc i 2 . By examining the behavior of the second mo- 
ment, we can clearly distinguish whether the system is 
in the diffusive regime or in the ballistic regime. To esti- 
mate the diffusion coefficient we discard the data in the 
ballistic regime. We also do not consider the time scale 
where the wave packet reaches to the edge of the whole 
system. Within the time scale considered in the present 
paper t < 200h/V, we do not observe any sign of the 
saturation of the second moment due to the finiteness of 
the system. To obtain the conductivity, we need the den- 
sity of states too. We estimate it by the Green function 
methodSi for strips the width of which is 12 < M < 30. 



III. NUMERICAL RESULT 

Let us first discuss how the diffusion of an electron in 
random magnetic fields depends on the strength of the 
fluctuation of random magnetic fields. We estimate the 
diffusion coefficients from the behavior of the second mo- 
ment shown in Fig. 1. It is clearly seen that the diffusion 
coefficient becomes smaller as the fluctuation of the ran- 
dom magnetic fields increases. The density of states at 
E = -0.5F is estimated to be p w 0.178, 0.181, 0.179, 
and 0.179(l/a 2 F) for h l{ = 0.2, 0.25, 0.3 and 0.4, re- 
spectively. It depends on the strength of the RMF very 
weakly near the band center. On the other hand, the 
diffusion coefficient D depends strongly on the strength 



of the RMF which is estimated to be AD = 85.22 ± 0.03, 
54.44 ± 0.04, 37.40 ± 0.03 and 21.30 ± 0.02(a 2 V/h) for 
h T{ = 0.2, 0.25, 0.3 and 0.4, respectively. In Fig. 2, the 
conductivity, evaluated from these diffusion coefficients 
and the density of states, is plotted for small values of 
h r { . Here we clearly see that the conductivity is inversely 
proportional to the square of h T t, namely, a oc \/h 2 f . 
This fact indicates that the present regime would be well 
described by the Born approximation. 15 In the recent 
experiment,- it is observed that the relative change of 
the resistivity by the random magnetic fields is in pro- 
portion to the square of the fluctuation of the random 
magnetic fields, which is qualitatively consistent with the 
present results (Fig. 2). 

Next, we apply an additional uniform magnetic field 
and examine its effect on the electron diffusion. The 
uniform magnetic field through each plaquette of the 
square lattice is denoted by </> U ni- The total magnetic 
field per plaquette is then </> un j + cj>, where un j is com- 
mon to all the plaquette and <f> is distributed indepen- 
dently as ©■ The second moments for h T i = 0.2 under 
various values of the uniform magnetic fields are shown 
in Fig. 3. Here it is clear that the diffusion of electrons 
is strongly suppressed by the additional uniform mag- 
netic field. Note that the strength of the uniform field 
considered is smaller than the fluctuation of the random 
field ^fWJWl = htf/VU = 0.0577 . . .. The diffusion 
coefficients and the density of states estimated for vari- 
ous values of the uniform fields are summarized in Table 
1. Here we see again that the density of states is not 
sensitive to the uniform magnetic field. We also perform 
numerical calculations for h T i = 0.3 and 0.4. In Fig. 4, 
the estimated conductivity is plotted as a function of the 
uniform field scaled by y/ (4> 2 ). For these values of the 
random magnetic field, it is commonly observed that the 
negative magnetoconductance occurs when the uniform 
field is weaker than the fluctuation of the random field. 
It is also observed that at 4> un i « W ((f> 2 ), the conductivity 
takes a value on the order ofe 2 /h and is, interestingly, in- 
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TABLE I: The diffusion coefficients D and the density of 
states p for h T f = 0.2 
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FIG. 3: Electron diffusion for ^ un i/^o = 0, 0.01, 0.02, and 
0.04, where </> un i stands for the strength of the uniform flux 
per plaquette. 



sensitive to the magnitude of the random magnetic field. 

In order to see whether this peak structure around zero 
field can be observed in the presence of the random po- 
tential (W 7^ 0), we also examine the electron diffusion 
with the additional random potential. We consider the 
case where h r { = 0.2 and W = 0, 1 and 2 and find that 
the peak structure disappears at W — 2 (Fig. 5). The 
random potential W — 2 on the square lattice yields 
the mean free path on the order of 4 lattice constants. 22 
This result indicates that to observe this enhancement 
of conductivity at zero uniform field we need fairly clean 
samples. 

Estimation of the mean free path in the random mag- 
netic field, especially for W — 0, is a subtle problem. 23 
If we define the relaxation time for the transport^ 3 , by 
Tt r = 2D/vp, where vp denotes the Fermi velocity, the 
mean free path Z r f in the random magnetic field can be 
defined by Z r f = vpT^ = 2D/vf- Near the band center, 
the average Fermi velocity can be estimated to be on the 
order of 2aV/h. The mean free path l T { is then evaluated 
to be about 21a, 9a, and 5a for h T i = 0.2, 0.3 and 0.4, 
respectively. 




FIG. 4: Magnetoconductance for htt = 0.2, 0.3 and 0.4. The 
uniform magnetic flux <f> un i is scaled by the strength of the 
random magnetic flux -J (4> 2 ). Inset: The normalized con- 
ductivity a /(To plotted as a function of </> un i ■ do. From the 
top, results for h T f = 0.4, 0.3 and 0.2 are presented. The 
normalized data do not fall on the one curve, suggesting the 
deviation from the Drude behavior. 
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FIG. 5: Conductivity in units of e 2 /h for h T f = 0.2 and 
W/V = 0, 1 and 2. 

IV. DISCUSSION 

We have demonstrated that in the absence of the uni- 
form magnetic field, the conductivity is inversely pro- 
portional to the square of the magnitude of the random 
field, and hence is very sensitive to the strength of the 
random magnetic field. In contrast, we have also shown 
that the conductivity takes a value on the order of e 2 /h 
and is insensitive to the strength of the random magnetic 
field if the uniform field is set to be 4> un [ fts y/{<fr 2 ) (Fig. 
4) . These two properties yield the peak structure at zero 
field, especially, in the case of the weak random magnetic 
fields. 

It may be useful to discuss our results in view of the 
Drude formula a = <?o/(l +lu^t 2 ), where u> c and r are 
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the cyclotron frequency and the relaxation time, respec- 
tively. Our results indicate that a oc r oc (0 2 ) _1 (Fig. 
2). The Drude theory then yields the form <j/<jq = 
1/(1 + A((j> un i ■ (Jo) 2 ), where A is a constant fixed by the 
density of electrons and independent of un i and 4>. We 
have found th at alt hough the data for small magnetic 
fields un i <C W (4> 2 ) seem to be consistent with this scal- 
ing form, it is unlikely that the Drude formula accounts 
for the behavior around un i w y/ {(j) 2 ). The deviation 
from the Drude theory, in which the resistivity is inde- 
pendent of the magnetic field, comes from the fact that 
the conductivity is insensitive to the magnitude of the 
random magnetic field at un j « yj (<j) 2 ) (Fig. 4). 

Let us consider this insensitivity of the conductiv- 
ity to the magnitude of the random magnetic field at 
</>uni ~ ■J {4> 2 ) m more detail. With this condition 
4>uni ~ yj (<fi 2 )i t ne system is almost equivalent to a sys- 
tem having the random magnetic fields distributed in the 
range < 4>/4>o < h and having no uniform magnetic field 
</>uni = 0. We have then evaluated the conductivity in 
such a system for 0.2, 0.3, 0.4 and 0.5 and, indeed, found 
that for all these values of h the conductivity is insensitive 
to the value of h and falls in the range 1.4 ~ 1.9(e 2 /h). 
This would be one of the significant transport properties 
specific to the random magnetic fields. 

In order to consider a possible relationship to the ex- 
periments, it may be useful to identify the correlation 
length of the random fields in experiments with the lat- 
tice constant of the present model. Our calculation 
then implies that the sample having mean free paths 
longer than 4 correlation lengths of the random fields 
is very sensitive to the application of uniform magnetic 



fields. This implication seems to be consistent with the 
experiments 16 where the dip structure is observed in sam- 
ples having a mean free path larger than the correlation 
length of the random magnetic field. 

In summary, we have investigated the electron diffu- 
sion in the random magnetic field with additional uni- 
form magnetic fields by the equation of motion method. 
We have found that a sharp peak at zero uniform field 
appears in magnetoconductance in the absence of the 
spatial correlation of random magnetic fields, where the 
semi-classical theory can not be applied. The width of 
the peak turns out to be on the order of the fluctua- 
tion of the random magnetic field. The conductivity at 
0uni ~ yj (4> 2 ) is found to be insensitive to the magni- 
tude of the random magnetic field. This peak structure 
disappears when the mean free path becomes shorter by 
introducing the random potential. Although the present 
method enables us to simulate very large systems, we can 
obtain the longitudinal conductivity a xx only. Detailed 
analysis of <j xx in experiments is required to examine the 
relevance of the present results. 
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